Theory of optical conductivity in detwinned YBa 2 Cu 3 
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We study theoretically the chain and in-plane optical conductivities of detwinned YBa2Cu3 !/ . 
We elaborate that the chain is superconducting for y > 6.67 while insulating for y < 6.67 due to the 
competition between the plane-chain coupling and the antiferromagnetic order, corresponding to a 
superconductor-insulator transition. Stemming also from the coupling between the plane and chain, 
a new peak emerges at a low frequency in the in-plane spectra in the superconducting state, while 
it disappears in the normal state. Our scenario accounts satisfactorily for very recent experiments. 

PACS numbers: 74.72.Bk, 74.25. Gz, 74.25.Jb 



YB^CnsOy (YBCO) is one of the most studied high- 
T c superconducting materials. A primary distinction be- 
tween this system and other cuprates is the presence 
of the quasi-one-dimensional CuO chain. Though it is 
widely believed that the main physics in this material is 
within the CuC>2 plane, the study of electronic properties 
in the chain and their influence on the superconducting 
plane may enable us to have a profound understanding 
on superconductivity in the present system 0, 0> 0, El- 
One of powerful tools to detect the electronic structure 
is the measurement of optical conductivity. It was re- 
ported that a pronounced a — b axis anisotropy exists 
in the spectra d, IE S H 0> \E1 E3 > which are strongly 
enhanced in the chain direction. By subtracting the a- 
axis spectra from the 6-axis (the chain direction) spectra, 
one can obtain experimentally the spectra contributed by 
the chain. In this way, the chain electronic structure re- 
vealed by the optical experiments shows some unusual 
features. Several experiments showed that the optical 
conductivity in the CuO chain exhibits a peak at a fi- 
nite frequency 0, El an d approaches to zero in the dc 
limit, implying that the chain is insulating though the 
Cu02 plane is still in the superconducting state. On 
the other hand, it was also reported that the chain spec- 
tra are characterized as a Drude-like peak located at the 
zero frequency, demonstrating the chain superconductiv- 

ity [111 US- 

Recently, a detail investigation of the infrared response 
in the detwinned YBCO [Tl| provided us a systematic 
doping dependence of the electrodynamics in the CuO 
chain. The optical spectra contributed by the chain dis- 
play a dominant Drude-like peak in the y = 6.75 sample. 
While as y decreases to 6.65, the Drude-like peak will 
shift to a narrow resonance peak at the finite frequency, 
indicating that the superconductor-insulator transition 
occurs in the CuO chain. In fact, the appearance of this 
transition around y = 6.65 may also be inferred from 
the doping dependence of the superfluid density 01 • A 
theoretical explanation of the chain superconductivity 
based on the plane-chain coupling was given by Morr and 
Balatsky [T^. however, an elaboration of the observed 



superconductor-insulator phase transition and the insu- 
lating nature of the chain is still awaited. In the mean- 
time, another important issue is how the in-plane spectra 
are influenced by the chain. Very recently, an experiment 
on detwinned ortho-II phase YBa2Cu306.5 observed an 
extra strong peak at 180 cm -1 (about 20 meV) in the 
a-axis optical conductivity in addition to the usually ob- 
served Drude-like peak and the mid-infrared (MIR) com- 
ponent in twinned samples |l3[ . Therefore, a consistent 
accounting for the optical response in both the plane and 
chain of the detwinned system is appealing. 

In this paper, we demonstrate that the competition 
of the plane-chain coupling and the antiferromagnetic 
(AF) order in the chain gives rise to the superconductor- 
insulator transition in the CuO chain, namely, the prox- 
imity superconductivity in the chain is induced by the 
plane-chain coupling when there is no AF order or it is 
negligible, while the insulating gap in the optical con- 
ductivity is caused by the AF order when its magnitude 
is appreciable. Starting from a self-consistent mean-field 
treatment for both the planar and chain Hamiltonian, we 
extend a simple existing model of the plane-chain 

coupling to include the emergency of the AF order in the 
chain at low dopings. Interestingly, we find a step-like 
rise of the AF order at the oxygen content y — y c = 6.67 
which is almost the same as the superconductor-insulator 
transition point observed experimentally |ll| . In the 
meantime, the chain optical conductivity shows an in- 
sulating gap below y c ; while above y Cl the chain optical 
conductivity rapidly acquires a Drude-like peak at the dc 
limit and shows a proximity-induced superconductivity. 
On the other hand, we find that a new peak occurs at 
a low frequency (0.18 J w 23mcV with J as the AF ex- 
change integral) in the in-plane optical conductivity due 
to the coupling between the plane and chain. These re- 
sults agree quantitatively with the recent experimental 
measurements 0, 0, and thus give a consistent pic- 
ture of the chain and in-plane optical conductivity in the 
detwinned YBCO. 

We start with a Hamiltonian which describes a system 
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with a plane and a chain per unit cell, 

H = H p + H c + Hj, (1) 
where H p describes the CuC>2 plane: 

h p = -*p E +h.c.-i E <%K+ h - c 
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(2) 



i/ c describes the y-direction CuO chain: 
ff c = -i c ]T C S S c CT + /i.c. + J c E(S 4 c .S J c -i<n J c ), (3) 

and iJj is the coupling between the plane and chain: 

Hj = -t^i^ + h.c). (4) 

Here (ij) denotes the nearest-neighbor (NN) bond and 
(ij)' the next NN bond. 

First, we use the slave-boson mean- held theory to de- 
couple the Hamiltonians (2) and (3). Then, the coupling 
of the planar fermions to spin fluctuations are included 
via the random phase approximation (RPA). In the slave- 
boson approach, the creation operators of the planar and 
chain electrons c^ c ^ are expressed by slave bosons 
carrying the charge and fermions ff£ representing the 
— ff} c ^b?^ c '. The mean-field order param- 
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spin, c, 

eters are defined as A? ■ 

(=b depend on if the bond (ij) is in the x or the y di- 

rection), X f ] = {f£° )] ffh = X P (c). Notice that we 
do not assume any superconducting pairing in the chain 
at the mean-field level. The AF order in the chain is 
defined as m c = (— l) l {Sf z ). At low temperatures we 
are concerned with, the boson condensation is assumed 

bf c) - (b? c) ) 
in the plane (chain). 

After Fourier transformation, the mean-field Hamilto- 
nian of the planar and chain fermions can be written as 



b";^' — * (ti- c ) = \/Sp(c), where <5 p ( c ) is the hole density 



+2N p j'( X 2 p + Al), 



(5) 



He — E( £ k/Jr/kcr + £ k+Q r /k+Q^o-/k+Q r cr) ~ 2J e ™ c 
kcr 

E ^(/kt/k+Qea + h - c -) + 2N c Mxl + ml), (6) 

ka 

where the summation over k in the chain Hamiltonian 
is in the magnetic Brillouin zone (MBZ) due to the AF 
order, i.e., —n/2 < k y < tt/2. Q c is the one-dimensional 
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FIG. 1: The hole density 8 C (dashed line) and AF order 
m c (solid line) in the chain versus the oxygen content y in 
YBa 2 Cu 3 !/ . 



AF wave vector with Q c = w/2. are given by e£ = 
—2(S p t p + J p x P )(cosk x + cos k y ) — 4<5 p £ p cos k x cos k v — 
fip and e£ = —2(5 c t c + J c x c ) cos k y — /i c , respectively. 
Correspondingly, the coupling of the planar and chain 
fermions can be written as 



Hi=t±J2(f^ + h.c). 



(7) 



k<r 



The renormalized Green's functions for the planar 
fermions G p due to the coupling to spin fluctuations are 
calculated by Dyson's equation in the Nambu represen- 
tation: 

G p (k, ilu)- 1 = G p0 (k, ilu)- 1 - S p (k, iw), (8) 
where the self-energy can be obtained from |l4j , 



G p0 (k - q, iu - iw m )a 3 , 



(9) 



G p o is the bare Green's function obtained from the mean- 
filed Hamiltonian and Xpi^aj^m - > w m +i6) is the RPA- 
type spin susceptibility 14]. The Green's function for 
the chain fermions G c in the Nambu representation is 
obtained directly from the above mean-field Hamiltonian. 

The real part of the optical conductivity a±(u>) is given 
by ciaa (w)p( c ) = -Imn aa (w) p { c )/u (a =X,y). Here the 
imaginary part of the current-current correction function 
Imn QCe (cL>) p ( c ) is expressed as 



[/(« + «') 



ImIJ aa (u) p(c) = EtT / ^'["aWftc)] 5 

k ly p(c) J 

-/(a/)]Tr[4 (c) (k, u + u')A p(c) (k, «')] , 

where v a is the a-component of the quasiparticlc 
group velocity and ^4 p ( c ) (k, to) is the spectral function 
[A p ( c )(k,w) = -(l/7r)ImGp( c) (k,cj)]. 

In order to investigate the doping dependence of the 
optical conductivity, we need to determine the planar and 
chain doping density 5 p ( C ) for a given oxygen content y 
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in YBa2Cu30 a . We determine 5 P versus y by usin g an 
empirical relation proposed recently by Liang et al. [l6j . 
On the other hand, in the parent compound (y = 6), the 
valence of Cu ion in the CuO chain is +1, and thus the 
electron density in the chain is 2. Correspondingly, we 
may have a conservation condition 2 — n c + 26 p = 2(y — 6), 
which relates the chain electron density n c and the pla- 
nar doping density S p with y. From this condition, we 
can derive the chain doping density S c = 1 
function of y, as shown in Fig.l. With the hole doping 
density, we can then solve the set of self-consistent equa- 
tions for mean-field parameters Xp(c), A p , m c , and n p ( c )- 
The input parameters are, t c = t. p = 2J, t p = — 0.45i p , 
J c = J p — 3 J/8, tj_ — 0.1 J. The AF coupling constant 
J in the plane is used as the energy unit (J ~ 130 meV). 
The chain AF order m c as a function of y at tempera- 
ture T « is shown as the solid line in Fig.l. We can see 
that the chain AF order m c emerges with the decrease of 
the chain doping density and reaches its maxima value 
at y — y rn — 6.61 where the chain doping density is 
zero(half- filled) 01 • We note that the AF order has a 
step-like rise from a negligible value to about a half of its 
maxima at y = y c = 6.67. 

The chain optical conductivity for different oxygen 
content y is plotted in Fig. 2. As y is increased to exceed 
y c = 6.67 [Fig. 2(a)], a prominent Drude-like peak at the 
dc limit shows up in the spectra, which indicates that the 
chain is superconducting. When the oxygen content y is 
decreased to y = 6.67, the Drude-like peak disappears 
and a finite frequency peak at 0.35J ~ 46meV occurs, 
leaving an excitation gap between this peak and the dc 
limit. This indicates that the system becomes insulat- 
ing at this doping density. With further decrease of the 
oxygen content, the finite frequency peak will move to a 
higher energy, leading to a larger excitation gap. Mean- 
while, the intensity of the spectra is suppressed heavily as 
shown in the inset of Fig.2(b). So, the superconductor- 
insulator transition occurs at y = 6.67, and the corre- 
sponding critical doping density in the chain is about 
gch ^ o.li. This result is consistent quantitatively with 
the recent experimental data |Tlj . 

From Fig.l, one can see that the AF order is negli- 
gible or disappears (for y > 6.9) in the doping range 
where a Drude-like peak in the optical spectra is observed 
and has an appreciable value when the optical spec- 
tra show an insulating behavior. More importantly, at 
y = y c = 6.67 where the insulator-superconductor transi- 
tion occurs, the magnitude of AF order has a step like in- 
crease, i.e., it rises from a negligible value to nearly a half 
of its maximum. This shows clearly that the transition is 
associated with the emergency of the AF order. The oc- 
currence of the AF order leads to a gap 2J c m c ks 24.4meV 
in the single-particle energy band of the chain, i.e., a 
48.8meV excitation gap in the particle-hole excitations. 
This excitation gap suppresses the proximity effect com- 
ing from the coupling to the superconducting plane which 
is otherwise effective for y > y c . Therefore, it is the com- 
petition between the AF order and the proximity effect 
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FIG. 2: The chain optical conductivity as a function of fre- 
quency in YBa2Cu30 H for different oxygen content y. The 
inset of Fig.(b) shows the spectra at y = y m = 6.61 where the 
calculated chain doping density is zero(half-filled). 

that leads to the superconductor-insulator transition in 
the chain. 

When m c equals to zero, the Hamiltonian for the sys- 
tem can be written as 4 x 4 matrix |14L fl5| . The gap 
symmetry induced by the proximity effect may be exam- 
ined by looking at the abnormal Green's functions of the 
chain electrons, F c (k,u>) — J^i ^3i(k, w)L74i(k, uj)/[iu> — 
E t (k,Lu)}, with Uij(k,uj) = A^? ± A jis ./F jk and A£ the 
superconducting gap in the plane with a d x 2_ y 2 symme- 
try |14| . Thus, the induced gap of the chain fermions is 
also of d x 2_ y 2 symmetry. 

Now we study how the in-plane optical conductivity is 
influenced by the coupling to the chain. The real part 
of the in-plane optical conductivity a\ (uj) is shown in 
Fig. 3 for the superconducting [Fig.3(a-b]) and the nor- 
mal state [Fig. 3(c)], respectively. In the superconducting 
state, a Drude-like peak at the dc limit and a MIR hump 
around oj ~ J is evident from Fig. 3 (a) and (b), which 
reproduces what has been observed in the twinned sam- 
ples 0. These features are also consistent with previ- 
ous theoretical calculations [l9l l20l f21 | in which only the 
in-plane electrodynamics is considered. A new feature 
observed here is that an extra peak emerges between the 
Drude-like peak and the MIR hump. Moreover, this extra 
peak disappears completely in the normal state, as seen 
from Fig. 3(c). These results are consistent with the very 
recent experimental data 0. We also note that there 
is only a weak a — b anisotropy in the spectra around 
this peak. In the meantime, the in-plane spectra in the 
superconducting state are nearly isotropic both in the dc 
limit and in high frequencies. As a result, the plane-chain 
coupling does not cause an obvious in-plane anisotropy in 
the optical conductivity. Therefore, we expect that the 
anisotropy observed in experiments, which is obtained by 
subtracting the a-axis spectra from the &-axis spectra, is 
mainly contributed by the chain contribution. 

Due to the coupling between the plane and chain, the 
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FIG. 3: Panels a)-c) plot the real parts of the in-plane op- 
tical conductivity as a function of frequency lu for different 
y. Panel a) shows the spectra for y = 6.75 sample in the su- 
perconducting state, panels b) and c) are the spectra for the 
ortho-II phase y — 6.5 sample with the planar doping density 
S p — 0.097 JLfJ in the superconducting state at T — 0.0005J 
and in the normal state at T = T c = 0.033J, respectively. 
Panel d) shows the quasiparticle energy E±k along the diag- 
onal direction in YBa2Cu3C>6.75. 



energy band of quasiparticle is split into two branches 
with frequencies E±k |14| . In the superconducting state, 
the low-energy optical response comes mainly from the 
charge excitations around the nodal direction due to the 
presence of the superconducting gap. In Fig.3(d), we 



plot E±k with k along the diagonal direction. Around 
the Fermi wave vector, there is a minimum u) m w 0.18 J 
in 2£+k, and the band -E_k is below this minimum. When 
the excitation frequency is below u> m , the optical response 
comes only from the -E-k- When ix> > Lu m , an additional 
scattering channel from the band -E+k contributes to the 
response and leads to a peak around Lu m . In the normal 
state, the excitations near the entire Fermi surface are 
available, so that the minimum in E + varies at different 
wave vectors. As a result, the combined effect of these 
excitations causes no extra peak. 

To summarize, we have elaborated that the competi- 
tion of the plane-chain coupling and the antiferromag- 
netic order in the chain leads to the superconductor- 
insulator transition in the optical conductivity of the 
CuO chain. Specifically, the proximity superconductiv- 
ity in the chain is induced by the plane-chain coupling 
while the insulating property is caused by the antiferro- 
magnetic order. Meanwhile, we find an extra peak in the 
in-plane optical spectra at low frequencies in the super- 
conducting state, which is between the Drude-like peak 
and mid- infrared component. This peak disappears in 
the normal state. These results are well consistent with 
the recent optical conductivity measurements. 
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